TWO-PARAMETER STOCHASTIC CALCULUS AND MALLIAVIN'S 
INTEGRATION-BY-PARTS FORMULA ON WIENER SPACE 



J.R. NORRIS 

Abstract. The integration-by-parts formula discovered by Malliavin for the Ito map on 
Wiener space is proved using the two-parameter stochastic calculus. It is also shown that 
the solution of a one-parameter stochastic differential equation driven by a two-parameter 
semimartingale is itself a two-parameter semimartingale. 



1. Introduction 

The stochastic calculus of variations was conceived by Malliavin [6j [8] as follows. 
Let (-Zf)t^o denote the Ornstein-Uhlenbeck process on Wiener space (W, W, fj) and let $ : 
W — > M. d denote the (almost-everywhere unique) Ito map obtained by solving a stochastic 
differential equation in M. d up to time 1. Then (zt)t>o is stationary and reversible, so, for 
functions /, g on R d , setting F = /oii) 1 G = go$ ) 

(1) E [{F(zt) - F(z )}{G(z t ) - G(z )}} = -2E [F(z ){G(z t ) - G(z )}] . 

Once certain terms of mean zero are subtracted, a differentiation of this identity with 
respect to t inside the expectation is possible, and leads to the integration-by-parts formula 
on Wiener space 

(2) / V i /($)r^'V i5 ($) C Z/x = - / fWLGdp, 

Jw Jw 

where LG and the covariance matrix T will be defined below. As is now well known, this 
formula and its generalizations hold the key to many deep results of stochastic analysis. 

Malliavin's proof of the integration-by-parts formula was based on a transfer principle, 
allowing some calculations for two-parameter random processes to be made using classical 
differential calculus. Stroock [TTJ 021 H3] and Shigekawa [TU] gave alternative derivations 
having a a more functional- analytic flavour. Bismut [1] gave another derivation based on 
the Cameron-Martin-Girsanov formula. Elliott and Kohlmann [3] and Elworthy and Li [5] 
found further elementary approaches to the formula. The alternative proofs are relatively 
straightforward. Nevertheless, we have found it interesting to go back to Malliavin's original 
approach in [8| and to review the calculations needed, especially since this can be done now 
in a more explicit way using the two-parameter stochastic calculus, as formulated in [9]. 

In Section [2] we review in greater detail the various mathematical objects mentioned 
above. Then, in Section El we review some points of two-parameter stochastic calculus from 
P|. Section [4] contains the main technical result of the paper, which is a regularity property 
for two-parameter stochastic differential equations. We consider equations in which some 
components are given by two-parameter integrals and others by one-parameter integrals. 
It is shown, under suitable hypotheses, that the components which are presented as one- 
parameter integrals are in fact two-parameter semimartingales. This is useful because 
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one can then compute martingale properties for both parameters by stochastic calculus. 
The sorts of differential equation to which this theory applies are just one way to realise 
continuous random processes indexed by the plane. See the survey [5] by Leandre for 
a wider discussion. But this regularity property makes our processes more tractable to 
analyse than some others. This is illustrated in Section [5j where we do the calculations 
needed to obtain the integration-by-parts formula. 

2. Integration-by-parts formula 

The Wiener space (W, W, //) over W 71 is a probability space with underlying set W = 
C([0, oo), M. m ), the set of continuous paths in W 71 . Let W° denote the cr-algebra on W 
generated by the family of coordinate functions w h- > w s : W — > M. m , s ^ 0, and let fi° 
be Wiener measure on W°, that is to say, the law of a Brownian motion in M. m starting 
from 0. Then (W, W, //) is the completion of the probability space (W, W°,/i°). Write W s 
for the ^-completion of a(w i— > w r : r ^ s). Let Xq,Xi, . . . ,X m be vector fields on M d , 
with bounded derivatives of all orders. Fix xq 6 M. d and consider the stochastic differential 
equation 

dx s = Xi(x s )dw l s +X (x s )ds. 

Here and below, the index i is summed from 1 to m, and d denotes the Stratonovich 
differential. There exists a map x : [0, oo) x W — > M. d with the following properties: 

• x is a continuous semimartingale on (W, W, (W s ) s ^o 5 A*)) 

• for ^/-almost all w € VF, for all s ^ we have 

x s (w)=x + / Xj(x r (u;))(9'u;* + / X (x r (t(;))(ir. 
Jo Jo 

The first integral in this equation is the Stratonovich stochastic integral. Moreover, for any 
other such map x' , we have x s (w) = x' s (w) for all s ^ 0, for yU-almost all w. We have chosen 
here a Stratonovich rather than an Ito formulation to be consistent with later sections, where 
we have made this choice in order to take advantage of the simpler calculations which the 
Stratonovich calculus allows. The Ito map referred to above is the map &(w) = x\(w). 

We can define on some complete probability space, (0, P) say, a two-parameter, con- 
tinuous, zero-mean Gaussian field (z s t : s,t ^ 0) with values in W n , and with covariances 
given by 

^{z? st z j s , tl ) = 5 l i{srKs')e-\ t - t '\l 2 . 

Such a field is called an Ornstein-Uhlenbeck sheet. Set zt = {z s t ■ s ^ 0). Then, for 
t > 0, both zq and zt are Brownian motions in M. m and (zo,zt) and (zt,zo) have the same 
distribution. We have now defined all the terms in, and have justified, the identity ([I]). 

Consider the following stochastic differential equation for an unknown process (U s : s ^ 
0) in the space of d x d matrices 

dU s = VXi(x s )U s dwi + VX (x s )U s ds, u = I. 

This equation may be solved, jointly with the equation for x, in exactly the same sense 
as the equation for x alone. Thus we obtain a map U : [0, oo) x W — > M. d ® (R rf )*, with 
properties analogous to those of x. Moreover, by solving an equation for the inverse, we can 
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see that U s (w) remains invertible for all s 0, for almost all w. Write £7* for the transpose 
matrix and set T s = U S C S U*, where 

C s = I U~ 1 Xi {x r )® U~ 1 Xi (x r ) dr. 
Jo 

Set also 

L s = -U s / U^Xiix^dwl + U s / U; l {V 2 Xi{x r )dw l r + V 2 X (x r )dr}T r , 
Jo Jo 

+ U s I U^WXiix^Xiix^dr 
Jo 

and define for G = g o $ 

LG = LiVtf(si) + Y'i'YX j >,[.<■ { ). 

We have now defined all the terms appearing in the integration-by-parts formula ([2]). We 
will give a proof in Section [5j 



3. Review of two-parameter stochastic calculus 

In [9], building on the fundamental works of Cairoli and Walsh [2] and Wong and Zakai 
|144 115], we gave an account of two-parameter stochastic calculus, suitable for the develop- 
ment of a general theory of two-parameter hyperbolic stochastic differential equations. We 
recall here, for the reader's convenience, the main features of this account. 

We take as our probability space (Q, J 7 , P) the canonical complete probability space of 
an m-dimensional Brownian sheet (w s t '■ s,t ^ 0), extended to a process (w s t ■ s,t £ R) by 
independent copies in the other three quadrants. Thus w s t = (u>l t , . . . , w™) is a continuous, 
zero-mean Gaussian process, with covariances given by 

E(ffi>^)=^(sAs')(iAt'), i,j = l,...,m, s,t^0, s',t'^0. 

It will be convenient to define also w;^ = si for all s,t € R. For s, t ^ 0, write T s t for 
the completion with respect to P of the cr-algebra generated by w ru for r G (— oo, s] and 
u G (— oo, i]. We say that a two-parameter process (x s t ■ s,t ^ 0) is adapted if x s t is 
J^t-measurable for all s,t ^ 0, and is continuous if (s,t) i— » x s t(u) is continuous on (M + ) 2 
for all oj € CI. The previsible cr-algebra on Cl x (M + ) 2 is that generated by sets of the form 
^4x (s, s'] x (t, t'] with A G T s t- If w e allow A G ^"soo in this definition, we get the s-previsible 
cr-algebra. 

The classical approach to defining stochastic integrals, by means of an isometry of Hilbert 
spaces, adapts in a straightforward way from one-dimensional times to two, allowing the 
construction of stochastic integrals with respect to certain two-parameter processes, in 
particular with respect to the Brownian sheet. Given an s-previsible process^ (a s (t) : 
s,t ^ 0), such that 

E / / a r (u) 2 drdu < oo 
Jo Jo 



We write any time parameter with respect to which a process is previsible, here s, as a subscript. Where 
previsibility is not assumed, here in t, we write the parameter in parentheses. 
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for all s,t ^ 0, we can define, for i = 1, . . . , m and all ti, ti ^ with ii i 2) one-parameter 
processes M and A by 



(3) 



Af, 



s rt 2 rs r-t 2 

I a r (t)d r d t w % rt , A s = / a r (t) 2 drdt. 

Jti JO Jti 



Then M is a continuous (^-" soo ) s ^o- mar tingale, with quadratic variation process [M] = 
A. A localization argument by adapted initial open sets (see below) allows an extension 
of the integral under weaker integrability conditions. By the Burkholder-Davis-Gundy 
inequalities, for all a £ [2, oo), there is a constant C(a) < oo such that 



(4) 



E 



S2 ftl 



Si Jtl 



a s (t)d s d t w 



St 



^ C(a)E 



S2 ft2 



Si Jt x 



a s (t) 2 dsdt 



a/2"* 



By an (s, t)-semimartingale, s-semimartingale, t-semimartingale, we mean, respectively, 
previsible processes {x s t : s,t^ 0), (p s t : s,t^ 0), (q s t : s,t ^ 0) for which we may write 



•Est %s0 



m 

E 

i=0 



XQt + XQO 

t 

I x~,. 



J 



m 

id r d u w l ru + } 

i,j=0 



1 Jo 



(x" u {r' ,u'))ijd rl d u w j r , u ) 



and 



Pst - Pot 



E 

i=0 ' 



(p' rt (v>'))id r d u >wl 



Qst — QsO 



E 



-1 Jo 



Wsu^idr'duWl'u- 



Here, (x' s ' t : s,t ^ 0) is a previsible process, having components (x' s ' t )j, subject to certain 
local integrability conditions, which are implied, in particular, by almost sure local bound- 
edness. The process (x" t (r,u) : s, t ^ 0, r, u 6 M) is required to be previsible in (u,s,t) 
and (Borel) measurable in (r, u), with x" t (r,u) = for r > s or u > i, and is subject 
to similar local integrability conditions. The inner and outer parts of the second integral 
are both cases of the stochastic integral at ([3]), or its t-analogue, or of the usual Lebesgue 
integral, and the value of the iterated integral is unchanged if we reverse the order in 
which the integrals are taken. The integrals appearing in the expression for x s t are called 
stochastic integrals of the first and second kind. The processes (p' st (u) : s,t ^ 0, u € M) 
and (q' st (r) ■ s,t ^ 0,r £ M) are required to be previsible in (u>,s,t) and measurable in u 
and r, respectively, with p' st (u) = for u > t and q' st (r) = for r > s, and are subject 
to similar local integrability conditions. For fixed t ^ 0, if (x s o : s ^ 0) is a continuous 
(•^so)s>o-semimartingale, then (x s t : s ^ 0) is a continuous (^ S i) s ^o-semimartingale, in the 
usual one-parameter sense. Also (p s t : s ^ 0) is a continuous (^ r s t) s ^o-semimartingale, for 
all t ^ 0. 
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The heuristic formulae 

m m ps pt 

d s d t x st = y2.{x" t )id s dtw\ t + V] / / (x" t (r,u))ijd s d u wl u d r d t w J rt , 



i=0 i,j=0 ' 

m rt 



pt 

dsPst = J2 (P'st(u))id s d u w l su , 
i=o 

dt?«t = J2 {q' s t(, r ))idrd t w l rt 



i=0 

provide a good intuition in representing the two-parameter increment 

d s d t X s t = X s+ ds,t+dt — %s,t+dt — X s+ ds,t + X st 

and the one-parameter increments d s p s t = p s +ds,t — Pst and d t q s t = Qs,t+dt — 1st in terms of 
a linear combinations of increments, and of products of increments of the Brownian sheet. 

By a (two-parameter) semi-martingale, we mean a process which is at the same time 
an (s, i)-semimartingale, an s-semimartingale and a t-semimartingale. Such processes are 
necessarily continuous. An (s, i)-semimartingale which is constant on the s-axis and t-axis 
is a semimartingale. By an obvious choice of integrands, the process (w s t : s,t ^ 0) is itself 
a semimartingale. The choice of lower limit —1 is useful to us in allowing as semimartingales 
a pair of independent M m -valued Brownian motions (z s q : s ^ 0) and (&ot : t ^ 0), given by 

ps pO pO pt 

z s o= / d r d u w ru , b ot = / d r d u w ru , 
Jo J -i J-i Jo 

which are moreover independent of (w s t ■ s,t ^ 0). Here and below, we bring one-parameter 
processes defined on the s or t axes into the class of two-parameter processes by extending 
them as constant in the second parameter. 

We say that a subset T> C (M + ) 2 is an initial open set if it is non-empty and is a union 
of rectangles of the form [0, s) x [0,t), where s,t ^ 0. A random subset PCSlx (IR + ) 2 is 
adapted if the event £ T>} is .^-measurable for all s,t ^ 0. For an adapted initial 

open set V, a process (x s t ■ (s, t) € I?) is a semimartingale in T> if there exists a sequence 
of adapted initial open sets V n j T>, almost surely, and a sequence of semimartingales 
(x™ t : s,t ^ 0), such that x st = x™ t for all (s,t) £ T> n for all n. The notion of an s- 
semimartingale in T> is defined analogously. We write C(^) f° r the boundary of T> as a 
subset of (1R+) 2 . In particular, if V = (M + ) 2 , then £(£>) = 0. 

The theory which we now describe is symmetrical in s and t. Where a statement is 
made for s, there is also a corresponding statement for t, which we shall often omit. Let 
(x s t : s, t ^ 0) and (x' st : s, t ^ 0) be s-semimartingales and let (a st : s, t ^ 0) be a 
locally bounded previsible process, for example, a continuous adapted process. There exist 
s-semimartingales which, for each t ^ 0, provide versions of the one-parameter stochastic 
integral and the one-parameter covariation process 

ps ps 
Cst = I a rtd r X r t, Cst = I d r x r td r x r 



lo Jo 



From now on, when we write these integrals, we assume that such a version has been chosen. 
We define also four types of two-parameter integral, each of which is a (two-parameter) 
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semimartingale. These are written 

pS pt pS pt 

Cst — / / &rudrduX rUl Cst — / / d r x ru d u y rui 

Jo Jo Jo Jo 

ps pt ps pt 

Cst — / / d r x ru d r d u y ru , Cst — / / d r d u x ru d r d u y ru . 
Jo Jo Jo Jo 

In the first and last integral, we require x to be an (s, t)-semimartingale, whereas, in the 

second and third, x should be an s-semimartingale. We require that y be a i-semimartingale 

in the second integral and an (s, t)-semimartingale in the third and fourth. All these 

integrals are defined as sums of certain integrals of the first and second kind with respect 

to the Brownian sheet. We refer to [9] for the details. We use the following differential 

notations: 



d s Zst = 


a s td s x s t 


means 




Zst - Zot = Cst) 


d s z s t = 


d s x std s x gt 


means 




Zst — Zot — L,st, 


d s dtz s t = 


a s td s dtx s t 


means 


Zst ■ 


~ ZgO — Zot + ^00 = Cst^ 


d s dtz s t = 


d s Xstdty s t 


means 


Zst ■ 


- Z S — Zot + Zoo — Cst) 


d s dtz st = 


d s Xstd s dtyst 


means 


Zst ■ 


- ZgO — Zot + Zoo = Cst, 


d s dtz s t = 


d s dtx s td s dty s t 


means 


Zst ■ 


~ ZgO — ZQt + Zqo = Cst- 


igrals Cst, 


C% and Cst an vanish if d s x s t = 


ctstds. It 


is shown in [9] that a 



identities hold among the various types of integral, which can be expressed conveniently 
in terms of this differential notation. Some identities assert the associativity of products 
involving a combination of three differentials or processes, the others are written as the 
following three rules 

d s (f(x s t)) = f{x st )d s x st + \f"{x s t)dsXstd s Xst, 

ds(a s tdtx s t) = dsdstdtXst + a s td s dtx s t + dsa s td s dtx st , 

d s (dtXstdty s t) = d s dtx s tdty s t + d t x s td s d t y s t + d s dtXstd s dty s t- 

These rules combine the usual calculus of partial differentials with Ito calculus in an obvious 
way. As a consequence, we can obtain a geometrically simpler Stratonovich-type calculus 
by defining, for processes (x s t '■ s,t ^ 0) and (y s t '■ s,t ^ 0), some further integrals, 
corresponding to the following differential rules 

X s td s X s t = XgtdYgt + ^d s Xstd s Yst, dgXgtdgYgt = dgXstdgYgt = dgXgtdgYst, 

where X s t may stand for any one of x s t,dtx s t,dtx s t and Y s t may stand for any one of 
y s t, d t y s t,d t yst- Then we have 

d s (f(x st )) = f'(xst)d s x s t, 
ds(a s tdtx st ) = d s a s tdtXst + a s td s dtXst, 
dg{d t Xstd t yst) = d s d t x s tdtyst + d t x s td s dty s t- 

The Brownian sheet {w s t '■ s,t ^ 0) and the boundary Brownian motions {z s o '■ s ^ 0) and 
{hot '■ t ^ 0) have some special properties, which are reflected in the following differential 
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formulae, for 1 ^ i,j ^ m, 

d s d t w l st d s d t w{ t = 5 lj dsdt, d s z l S Qd s z? sQ = 5 lj ds, dtb l ot d t bl t = 5 lj dt, 
and, for any semimartingale (x s t : s,t ^ 0), 

d s x st d s dtw l st = d t x s td s d t w % st = 0. 



4. A REGULARITY RESULT FOR TWO-PARAMETER STOCHASTIC DIFFERENTIAL 

EQUATIONS 

We discussed in [9] a class of two-parameter hyperbolic stochastic differential equations, 
in which there is given, for a system of processes (x st ,p s t,q s t '■ s,t ^ 0), one equation for 
the mixed second-order differential d s dtx s t, together with two further equations for the 
one-parameter differentials d s p s t and dtqst- We review briefly the details below, and then 
give a new regularity result, which we need for our application to Malliavin's integration- 
by-parts formula, but which may be of independent interest. This result concerns the 
process (p s t : s,t ^ 0) (and analogously also (q st : s,t ^ 0)), which, since integrated 
in s, has naturally the regularity of an s-semimartingale. The point at issue is whether 
{p s t : s, t ^ 0) is a full (two-parameter) semimartingale. A method to establish this is 
stated in [9l pp. 299, 315-316], but the argument given is incomplete. A full proof is given 
below in Theorem 14.21 As an illustrative example, we note that, if (w s t : s, t 0) is a 
Brownian sheet with values in W 71 , then the result will show that there is a two-parameter 
semimartingale (x s t : s, t ^ 0) such that, for all t ^ 0, the process (x s t : s ^ 0) satisfies the 
one-parameter stochastic differential equation 

d s x st = Xi(x st )d s w l st + X (x st )ds, 

with given initial values xot = xo, say. This is useful because, now, despite the irregular 
dependence of the Brownian sheet on t, we can use a differential calculus in t as well as in 
s. 

Consider the class of hyperbolic stochastic differential equations in (M + ) 2 of the form 

(5) d s d t x st = a(d s d t w st ) + b(d s x st , d t x st ), 

(6) d s p s t = c(d s x st ), 

(7) d t q s t = e(d t x st ). 

Here w s t = (wl t , . . . , w™),with (w\ t : s, t ^ 0), i = 1, . . . , m, independent Brownian sheets, 
as above. The unknown processes (x st : s,t ^ 0), (p s t '■ s,t ^ 0) and (q s t : s,t ^ 0) take 
values in M. d , M. n and M. n , respectively, and are subject to given boundary values (x s q : s ^ 0), 
(xo* : t ^ 0), both assumed to be semimartingales, and (pot ■ t ^ 0), (q s o : s ^ 0), both 
assumed continuous and adapted. The coefficients a, b, c, e are allowed to have a locally 
Lipschitz dependence on the unknown processes, with the restriction that b depends only 
on x. Thus, for example, we would write a(x s t,p s t, q s t, d s dtw s t) and b(x s t, d s x s t,dtx s t), but 
have not done so in order to keep the notation compact. Moreover, we allow a dependence 
on the differentials which is a sum of linear and quadratic terms. Thus, in an expanded 
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notation, we would write 

d s dtx st = a 1 (d s d t w st ) + a 2 {d s d t w st ,d s d t w st ) 

+ b u (d s x st , dtXst) + bi2(d s x s t, d t x st , d t x st ), 
+ b 21 (d s x st ,d s x st ,dtx st ) + b 22 (d s x st ,d s x st ,d t x st ,d t x st ) 
d s Pst = ci{d s x st ) + c 2 (d s x st , d s x st ), 
dtQst = ei{d t x st ) + e 2 {d t x st ,d t x st ), 
where, for k = 1, 2, 

ad iran . . iren , Trad * 



6 ife :R d ^R d ® ((R d )*)® J ' +fe , 

Cj - : R d x R" x R" -» R n ® ((]R d )*)^', 

e fc : R d x R n x R n -» R n <g> ((IR d )*f fc . 

We may and do assume with loss that a 2 ,b\ 2 ,b 2 i,b 22 ,c 2 ,e 2 are symmetric in any pair of 
repeated differential arguments. 

By a local solution of dSHZ]) with domain T> we mean an adapted initial open set T>, 
together with a semimartingale (x s t : (s,t) G £>), an s-semimartingale (p s t : (s,t) G 2?), and 
a i-semimartingale : (s,t) G T>), all continuous on V, such that, for all (s,t) G P, 

ps pt ps ft 

x s t = x s0 + x ot -x 00 + / a(d r d u w ru ) + / / b(d r x ru , d u x ru ), 

Jo Jo Jo Jo 



Pst =Pot+ / c(d r x rt ), 
Jo 

Qst = q s o + / e(d u x su ) 
Jo 



Given such a solution, for each t ^ 0, we can define processes (u s t '■ (s, t) G T>) and 
(u* t : (s, t) G V), taking values in R d x (R d )* and R d x x (R d )* respectively, by solving 
the linear one-parameter stochastic differential equations 

(8) d s u st = b u (d s Xst,-)ust + b 12 (dsX s t,d s Xst,-)ust, 
d s u* t = u~t{b 12 (d s x st ,u s r,u S f) 

(9) + b 22 (d s x s t , d s x s t ,u st -,u S f) - b n (d s x s t , b 12 (d s x s t , u s t ■ , u st •)) }■ 

Here denotes the inverse of the linear map u s t- For fixed t ^ 0, almost surely, u s t 
remains in the set of invertible maps while (s,t) G T>. To see this, one can obtain formally 
a linear equation for the process (v,g t : (s,t) G T>), and then check that its solution is 
indeed an inverse for u s t- Similarly, for each s ^ 0, we can define processes (v s t '■ (s,t) G T>) 
and (v* t : (s,t) G V), taking values in R d x (R d )* and R d x (R d )* x (R d )* , by solving the 
analogous equations 

(10) d t v st = bii(-,d t Xst)v s t + b 2 i(-,d t x sU d t x st )v st . 

dtv* st = vJt{b 21 (v s r,v s r,dtx s t) 

(11) + b 22 (v S f,Vsf,dtXst,d t Xst) - bii{b 2 i{vsr,Vsr,d t Xst),d t Xst)}- 



s 



We specify initial conditions u q = vqq = I, so determining completely (uq s : s ^ 0) and 
(t>ot : t ^ 0). Then we complete the determination of the above processes by specifying that 
u Qt = v ot , ul t = 0, v s0 = u s0 , and v* = for all s, t ^ 0. Let us say that (x st ,p st , q st : (s, t) £ 
D) is a regular local solutioifl if there exist continuous s-semimartingales (u s t ■ (s,t) G T>) 
and (u* t : (s,t) € V) satisfying, for each t ^ 0, the equations (JHH9|), and if there exist also 
continuous i-semimartingales (v s t : (s,t) € £>) and (t>* t : (s,t) € I?) satisfying, for each 
s ^ 0, the equations (jlOHllh . A local solution is maximal if it is not the restriction of any 
local solution with larger domain. The notion of a maximal regular local solution is defined 
analogously. We assume that the boundary semimartingales (x s q : s ^ 0), (xQ t : t ^ 0), 
(pot '■ t ^ 0) and (g s o : s ^ 0) are regwZajQ. By this we mean that the Lebesgue-Stieltjes 
measures defined by their quadratic variation processes and by the total variation processes 
of their finite variation parts are all dominated by Kds, or Kdt as appropriate, for some 
constant K < oo. We give a result first for the case where b = 0. 

Lemma 4.1. Assume that 6 = 0. Let U be an open subset of M. d x M. n x M n and let 
m : U — > [0, oo) be a continuous function with m(x,p,q) — > oo as (x,p,q) — > <9£/. Assume 
that, for all M ^ 1, i/te coefficients a,c,e are bounded and Lipschitz on the set Um = 
{(x,p, q) € J7 : m(x,p,q) < M}. Then, for any set of regular boundary semimartingales 
(x s0 : s ^ 0), (aj t : t ^ 0), (p t : * > 0) and (a s0 : s ^ 0), with (x o,Poo, <7oo) G ^, 
equations ([SHU) aaue a unique maximal local solution (x s t,Pst,Qst '■ ( s ^t) £ with values 
in U. Moreover, we have, almost surel 1 ^ 

sup m(x ru ,p ru ,q ru ) -» oo as (s, i)t 

Proof. In the case where m is bounded (so Um = U = M. d x IR" - x M n for large M), the 
existence of a (global) solution is proved in [9l Theorem 3.2.2]. The proof is of a standard 
type, using Picard iteration, GronwalFs lemma and Kolmogorov's continuity criterion, and 
gives also the uniqueness of local solutions on the intersections of their domains. When m 
is unbounded, we can find, for each M ^ 1, bounded Lipschitz coefficients om,cm,gm on 
M. d x W 1 x 1", which agree with a, c, e on Um- F° r each Mo #J 1, the corresponding global 
solutions (Xg{,Pst,qst : s,t ^ 0) agree, for all integers M ^ Mo, almost surely, on T>m , 
where 

P A/ = {M) G (M + ) 2 : sup rn(xfi,pM,qM)^M}. 

Hence, we obtain a local solution with all the claimed properties by setting T> = Um^m 
and by setting, for all M ^ 1, (x st ,p st ,q st ) = (x^,p^,q^) for all (s,t) € V M \ V M -l- □ 

Our main result deals with the case when b is non-zero. 

Theorem 4.2. Assume that the coefficients a, b, c, e are uniformly bounded and Lipschitz. 
Then, for each set of regular semimartingale boundary values (x s q : s ^ 0), {x$t : i ^ 0), 



It is not hard to see that, for any local solution, the processes just defined have previsible versions, 
which are then s-semimartingales or f-semimartingales, depending on the variable of integration. However, 
we have not determined whether they have a continuous version in general. 

No connection with the notion of regular local solution is intended. 

To clarify, we mean that, for all (s*,t*) G C,{T>), the given limit holds whenever (s,t) f (s*,t*). In 
particular, in the case where T> = (R + ) 2 , there are no such points (s*,t*) and nothing is claimed. 
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(Pot '■ t ^ 0), (q s o : s ^ 0), the system of equations dSHZ]) has a unique maximal regular 
solution, with domain T> say. As (s,t) f C(^)> we have 

(12) m si = sup \(u slt/ ,u~],,v s/t/ ,v~l)\ -> oo. 

Moreover, if c has Lipschitz first and second derivatives and has no dependence on q, then 
(p s t : s,t ET>) is a semimartingale in V . 

Proof. We consider first the question of existence. We follow, to begin, the strategy used in 
the proof of [HI Theorem 3.2.3]. Consider the following system of differential equations, for 
unknown processes y s t, z s t, x' st , u s t,u* t ,p s t, x" t: v s t,v* t , q s t, taking values in W d , M. d , W d , M. d <S> 
(R d )*,R d ® (R d )* <g> (R d )*,R n ,M ,i ) M d ® (R d )*,R d ® (R d )* ® (M d )*,M n respectively: 

(13) d s d t y st = u~^a(d s d t w st ) - u* st {u~^a(d s d t w st ) ® u~^a{d s d t w st )), 

(14) d s d t z st = v~ t l a(d s d t w st ) - vl t (v~^-a(d s dtw st ) ® v~^a{d s d t w s t)), 

(15) d s x' st = v st (d s z st + v* t d s z st ® d s 2r st ), 

(16) 4u st = &n(u st (d s z si + v* st d s z st <g> d s Zgt), O^st + b2i(v st d s z s t,v st d s z st , -)u st , 
d s u* st = u~t{bi2{v st {d s z st + v* st d s z st <g) d s z st ),u S f,u S f) 

(17) + b 2 2(v s td s z s t,v st d s z s t,u S f,u s f) - bu(v st d s z st , bi 2 {v st d s z st ,u S f,u S f))}, 

(18) d s p st = c(v st {d s z st + v* st d s z st <S> d s z st )), 

(19) d t x" t = u st (d t y s t + u* st d t y st ® dtUst), 

(20) dtv 3t = &u(-, u st (d t y s t + u* t d t y st <g> d t y s t))v s t + &12G, u st d t y s t, u st dty s t)v s t, 
dtv* st = v~ t l {b 2 i(v S f,v S f,u s t(d t y s t + u* t d t y st ® dty s t)) 

(21) + b 2 2(vsf,v S f,u st dty s t,u st dty s t) - bn(b 2 i(v S f,v S f,u s tdtyst),Ustdtyst)}, 

(22) dtg st = e(u st {d t y s t + u* st d t y s t <8> ckyst)). 

We evaluate the coefficients a, 6, c and e here at {x' st ,p s t,q s t) (rather than at a;^). Note 
that this system has the same form as the system ([5HZ1) with 6 = 0. We use the boundary 
conditions given above for u s t,p s t,v s t,q s t- Define boundary values for y s t and by 

(23) eZ s y s o = = ^c^o, d t y ot = d t z ot = u^d t x ot , Voo = ^oo = 0. 

Set Uq 4 = v* = and use the given boundary values (xpt : i ^ 0) for x' st and (x s o : s ^ 0) 
for x" t . Define, on the set U where u and v are invertible, 

m(y,z,x',u,u*,p,x",v,v*,q) = |(tt, u~\ v, u _1 )| + |(tt*,«*)|. 



Then the preceding lemma applies, to show that (|13H22[) has a unique maximal local solution 
with the given boundary values, with domain T> say, such that u st and v s t are invertible for 
all (s,t) € V, and such that, almost surely, as t | CC^Oj either 

(24) m st = sup |(« s 'f jitj,],,^ ,u^)| t 00, 
or 

(25) n st = sup \(u*, t ,, v*, t ,)\ I 00. 

s%s,t%t 
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Now v s t and v* t are continuous i-semimartingales (in T>) and z s t is a semimartingale. 
Moreover dta s td s dtz s t = for any t-semimartingale a s t- Hence, by [U Theorem 2.3.1], x' st 
is a semimartingale and we may take the t-differential in (|15p to obtain 

d s d t x' st = d t v st (d s z st + v* t d s z st <g) 

+ v st (d s d t z st + d t v* t d s z st (8) d s <Zst + v* t d s d t z st d s d t z st ) + d t v s t(d t v* t d s z st ® d s £ st ) 

= a(d s d t w st ) + b{d s x' st ,dtx" st ). 

Similarly, by taking the s-differential in (|19p . we obtain 

44#"t = a(d s d t w st ) + b(d s x' st , dtx" t ). 

We also have Xq = Xq and 

4x' s0 = v s0 d s z s0 = d s x" , d t x' ot = uotd t y t = d t XQ t , 

so x' st = x" st for all (s, t) G X>, almost surely. Denote the common value of these processes 
by x S f. Then (x s t : (s,t) G T>) satisfies ©. On using (USD and CED to substitute!! for 
and ci t y si in ([TBI [M [20l [22]), we see also that p st , q st , u st , u* st , v st , v* t satisfy ([6HII]) 
respectively. Hence (x s t,p s t, 1st '■ (s,t) G X>) is a regular local solution to (|5H7|), which is 
moreover maximal by virtue of (|24H25p . 

We turn to the question of uniqueness. Suppose that (x s t,Pst,<ist ■ G V) is any reg- 
ular local solution to ((SH7|)- Write {u s t,u* t ,v s t-,v* st : (s,t) G V) for the associated processes, 
satisfying ([H HTT]) . Define semimartingales (y st : (s,t) G V) and {z s t ■ (s,t) G T>) by 

(26) d s dtjjst = u~ta(d s d t w st ) - u* st (u~^a(d s d t w st ) ® uj t l a(d s d t w st )), 

(27) d s d t z st = v~ t l a(d s d t w st ) - v* st (v~^-a(d s dtw st ) v~ t l a(d s d t w st )), 

with boundary values (|23|) . The following equations may be verified by checking that the 
initial values and differentials of left and right hand sides agree 

(28) d s x s t = v st (d s z st + v* t d s z st ® d s z st ), d t x st = u s t(d t y st + u* t d t y st ® d t y st ). 

Then, using these equations to substitute for d s x s t and dtx s t in ([6 HTT]) . we see that 
(y s t,z s t,x s t,u s t,u* t ,p s t,x s t,v s t,v*j,q st : (s,t) G V) is a local solution to (H3H22D- By lo- 
cal uniqueness for this system, D C D and (x s t,Pst,q~st) = fast, Pst, Vst) for all (s,t) G V, 
almost surely. Thus (x s t,Pst,Qst '■ (s,t) G T>) is the unique maximal regular local solution 

to dHHZD- 

Our next goal is to obtain ath-moment and L a -H61der estimates on the process 
(x s t,p s t,q s t,Ust,ul t ,v s t,v* t : (s,t) G T>), for a G [2, oo). Write K for a uniform bound 
on a,b,c,e which is also a Lipschitz constant for b. Fix M, N, T 1 and set 

T>m = {(s, t) G V : s, t < T and m st < M}, 
T>m,n = {{s, t) G V : s,t ^ T, ?n st ^ M and n s t ^ iV}. 

Fix a and define 

g(s,t)= sup E(|(n*, t ,,<, t ,)| a l {{s / it / )g7?Miv} ). 



Such substitutions result in differential formulae corresponding to valid identities between processes. 
This is because the two-parameter stochastic differential calculus is associative, as mentioned above, and as 
discussed in [9] pp. 290-291]. 
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Let (a s : s ^ 0) be a locally bounded, (J r soo ) s j>o"P rev i s ibl e process. The following identities 
follow from equations (f2"7|) and (j28j) : for (s, i) G X>, respectively in and R rf <g) M d , 
(29) 



a r d r x r Q + 



o Jo 



a r v rt {v r ^a(d r d u w ru ) + (v* t - v* ru )(v r ^a(d r d u w ru ))® 2 } 



and 
(30) 



a r d r x rt ® d r x rt = \ a r d r x r0 <g> d r x r0 + / / a r (v rt v r ^a{d r d u w ru ))® 2 . 



o Jo 



Hence, using the estimate we obtain a constant C = C(a, K, M, T) < oo such that, for 
all s, t ^ 0, 



E 



(31) 

and 
(32) 



Qij> dtp 

€ CE 



1 



{{s,t)eD MtN } 

1/2 



o J o 



E 



a r cL 



l a r |(|^ri I + |v*J)<ir<iu 



| a r | dr 



1 



{(s,t)€X> M ,iv} 



1 



{{s,t)£D M>N } 



Here and below, we suppress any dependence of constants on the dimensions d, n, m. If we 
allow C to depend also on N, then (|3ip may be simplified to 



(33) 



E 



CLy df X'p-fc 



a 2 .dr 



a/2 



1 



{{s,t)ev MiN } 



We use these estimates, along with analogous estimates for integrals dtx s t, in the equations 
(|9|) and (fTTj) . to arrive at the inequality 



g(s, t) ^ C (l + J g(s', t)ds' + J g(s, t')dt' 



for a constant C = C(a,K,M,T) < oo. Since N < oo, we know that g(s,t) < oo for all 
s,t, so this inequality implies that g(s,t) ^ C for another constant C < oo of the same 
dependence. Similar arguments yield a further constant C < oo of the same dependence 
such that, for all s, s' ^ and all t, t' 0, 

(34) E(\(x st ,u st ,u* st ,p st ) - (x s > t ,u s/t ,u* s , t ,p s < t )\ a l{( St t),(s',t)£V M , N }) «S C|s - s'| q/2 
and 

(35) E(\(x st ,v st ,v* t ,q st ) - {x st ',v st ,, v* t ,, q s t')\ a l{( s ,t),(s,t')eV MiN }) < C|t-t'| a/2 . 
Here, we have used Cauchy-Schwarz to obtain in an intermediate step 

KJdrdu < |«- s'l 1 / 2 ( / / KJ^rdu^ . 



s ^0 
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On going back to 
constant C < oo of the same dependence such that 



(36) 
and 

(37) 
Now 



E 



E 



a r (d r 



i 



t),(s,t')ev M , N } ) ^C\t- t'\ a/2 I E 



ain, 


using (UD aj 


j;ain 








1/2 


(e 


/ a^ds 


1 






Jo 







i 



{(s,t),(s,t')eT> M ,N} 



(e 


/ a^ds 




Jo 



a\ 1/2 



d s (u s tu st >) = u st 1 {b(x ai /,d a x s t>,-) - b(x st ,d s x st ,-)}u st > 



bu(x s t,d s x st ,-)}u s t'- 



We have made explicit the dependence of b and b\\ on x st or x st >- We use the estimates 
Plj) . pZ]) . dHSHSZl) to find a constant C = C(q, M, T) < oo such that 



(38) 



E(K* - u st ,\ a l {{sMstt , )€VMN} ) < C|t - t'\ a ' 2 . 



Moreover, the same estimates, applied to the difference of ([9]) at t and at t', show that C 
may be chosen such that 



(39) 



- ^ri{( s ,t),( s ,t')e^,iv}) < C\t - t 



not/2 



Since C does not depend on N, by monotone convergence, we can replace T>m,n by Pm in 
these estimates By symmetry, there are analogous estimates for v s t and v* t . Hence, using 
[U Theorem 3.2.1], almost surely, for all M ^ 1, n s t remains bounded on T>m- Thus (|25p 
implies (|24|) so, in any case, (|12p holds. 

It remains to consider the case where c has Lipschitz first and second derivatives and has 
no dependence on q, and to show then that (p s t ■ (s,t) € T>) is a semimartingale. For ease 
of writing, we shall assume that c has no dependence on x either. This is done without 
loss of generality, by the device of adding to our system the equation d s x st = d s x st , thus 
making x s t a component of p s t- 

We seek to find a solution in a smaller class of processes, in which p s t is a semimartingale. 
Recall that 

(40) d s p s t = c(d s x st ) = ci{p st ){d s x st ) + c 2 (p s t)(d s x st , d s x st ). 

By Ito's formula, if p s t is a semimartingale, then 

d s d t pst = c'(d t p s t,d s x st ) + \c'{dtPst, d t p s t,d s x st ) + c(d s d t x st ) + c'(d t p s t, d s d t x st ) 
+ 2c 2 (d s x st , d s d t x st ) + 2c' 2 (d t p st , d s x st , d s d t x st ) 
= c'(d t pst, d s Xst) + lc"(dtp s t, d t pst,d s x st ) + c(a(d s d t w st )) + c(b(d s x st , d t x st )) 
+ c'(d t p st , b(d s x st , d t x st )) + 2c 2 (d s x st , b(d s x st , d t x st )) 
+ 2c 2 (d t pst, d s x st , b(d s x st , d t x st )). 
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Here we are writing c',c" for the derivatives with respect to p. We set d = d + n and 
combine this equation with the equation ([5]) to obtain a two-parameter equation for the 

M^-valued process x s % = ( Xst ) , which we can write in the form 



K Pst / 

(41) d s d t x st = a(d s d t w st ) + b(d s x st , d t x st ). 

(The ~ notation in this paragraph has nothing to do with that used in the paragraph 
on uniqueness above.) We impose regular semimartingale initial values x s q = (^p S °^j an< ^ 

XQ t = [ Xot ] , where (p s o : s ^ 0) is obtained by solving the one-parameter equation (HU1) 

along x s q. Introduce the two companion equations for d X d matrix- valued processes u s t 
and v s t 

(42) d s u st = bn(d s x st , -)u st + b 12 (d s x st , d s x st , -)u s t, 

(43) dtVst = bn(-,dtx s t)v s t + b 2 i(-,dtx s t,dtx s t)v s t. 

Impose boundary conditions for u st and v s t analogous to those for u s t and v st . Write ([7]) 
in the form 

(44) d t q st = e(d t x st ). 

By assumption, there exists a K' < oo which is both a uniform bound for a, b, c, e and 
is also a Lipschitz constant for b,c,c',c". We can then find a uniform bound K < oo on 
a, b, e, which is also a Lipschitz constant for b, and which depends only on K' . The above 
argument shows that the system of equations (|41H44p has a unique maximal regular solution 
(x st ,q s t,u s t,v st : (s,t) € V), with the property that, as (s,t) f Q{T>)-, almost surely, 



m st : = sup \(u ru , u ru , v ru , v ru 

r<s,u<t 



^)lt°C. 



Write 



and use analogous block notation for the tensors u* t and Note that 
where 

f(d s xl t ) = c{b{d s x\ t , ■)) + 2c 2 (d s xlt, b(d s x\ u ■)), 

g{d t x st ) = c'(d t x 2 t ,-) + \d'(d t x 2 st ,d t x 2 st ,-) + c(b(-, d t x\ t )) + c'(d t x 2 st , &(■, d*a4)). 

Here, we have written b(d s x s t,-) as a short form of &n (d s x s t, •) + 612 (d s x s t, d s x s t, •), and 
analogously for b(-,d t x s t) and b(d s x s t,-). On multiplying out in blocks, we see that the 
process (x* t , x^, q^, m^ 1 , (ii* t ) m , i^ 1 , (^) m : € £>) satisfies equations (JoT fTTj) . Hence, 
we must have and (a^ t , x 2 t , q st , u£, vlj) = (x st ,Pst,Qst,Ust,v st ) for all (s,t) 6 P. In 

particular, (p st : (s,t) € £>) is a semimartingale. 
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It remains to show that T> = T>, which we can do by showing that, almost surely, rh s t 
remains bounded on T>m,n — 'DC\'Dm,n, f° r an M, N ^ 1. We first obtain a Holder estimate 
in t for p s t. We have 

dsipst ~ Pst') = c(p st , d s Xst) - c(p st r,d s x st >), 
where we have now made the dependence of c on p explicit. Set 

f(s) = E (\ Pst - p s t'\ a l {{smS}tl)e -3 M>N} ) ■ 

We use the estimates (|32p and f)33[) to obtain a constant C = C(a, K' , M, N,T) < oo such 
that 

/(sK C (V - i'r /2 + f /(r)drV 

This implies that /(s) ^ C|£ — t'| a / 2 for all s ^ for a constant C < oo of the same 
dependence. We now know that, for such a constant C < oo, we have 

(45) e (W - P*r i {(s , t ),(^)«w) < c (i s - s 'i Q/2 + I' " ^ 

We turn to u s t and S S £. The following equations hold 

d s u]t = b(d s x st , ■)«]£, d t v l £ = b(-,d t x st )vl 2 , d t v 22 = g(d t x st )v 1 ^. 

By uniqueness of solutions, we obtain u^ 2 = u stUQ^u^ so, in particular, itj, 2 . = 0. Similarly, 
vl$ = v s tv~QvlQ, so = 0. Since {tot = vot and u s q = v s q, we deduce that u\ 2 = v\ 2 = 0. 
Then dtv 22 = 0, so v 22 = v 2 q = u 2 ^. We also have the equations 

d s v%t = f(d s x s t)u st + c'(., d s x st )u 2 J, d s u 22 = c'(., d s x st )u 22 , d t v 2 \ = g(d t x s t)v st 
and we note that 

— i _ f «rt ^ —i-f v 7t o A 

" \-{uf t )- i uf t u- t i (u 22 )-i) ' ^ - U^r^v* 1 k 2 ? rv ' 

and 

d s(«s 2 ) _1 = -(« 22 )~ 1 C / (.,4^) + (uf)" 1 ^.,^^)^.,^^). 

We use the inequalities (|32p . (|33|) and (|45p . and an easy variation of the argument leading 
to and (J3"8j) to obtain a constant C = C(a, K', M, N, T) < oo such that 

(46) e (Ktw.tvj) - (fi rt ,fi^)r itotwy.iorf^}) < c?(i« - ,r /2 + 1* - tr /2 )- 

Then, using [91 Theorem 3.2.1] as above, we can conclude that, almost surely, (u s t,u~^~) 
remains bounded on T>m,n- It remains to show that the same is true for (vstjV^ 1 ) and, 
given the relations already noted, it will suffice to show this for v 2 l . We have 

d s u* st = u~t{bi 2 {d s x s t,u s t-,u st -) 

+ b 2 2(dsXst,d s x s t,u S f,u S f) - bii(d s x s t,bi2(d s x st ,Usf,u S f))} 

= h(x st ,p s t,u s t,u~ t ,d a x st ), 

where h is defined by the final equality and where we have used © to write d s x s t in terms 
of d s x s t- A variation of the argument used for u st shows that, almost surely, u* t remains 
bounded on T>m,n- Then, we can use the ~ and t-analogue of equations (I29p and (130)) to 
express v 2 } as a sum of integrals with respect to (xot)Pot ■ t ^ 0) and (w s t : s,t^ 0). This 
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leads, as above, to L a -Holder estimates which allow us to conclude that, almost surely, v^j 
remains bounded on T>m,N, as required. □ 

5. Derivation of the formula 

Let (w s t '■ s,t 0) be an M m -valued Brownian sheet and let (z s q : s 0) be an inde- 
pendent IR m -valued Brownian motion. Thus w s t = (wl t , . . . , w^) and z s q = (z] , . . . , z^), 
and each component process is an independent scalar Brownian sheet, or Brownian motion, 
respectively. The two-parameter hyperbolic stochastic differential equation 

(47) d s d t z st = d s d t w st -\d s z st dt, s,t^0, 

with given boundary values (z s q : s ^ 0) and Zot = 0, for t ^ 0, has a unique solution 
(z s t : s,t ^ 0). Set z t = (z s t : s ^ 0), then (zt)t^o is a realization of the Ornstein-Uhlenbeck 
process on the m-dimensional Wiener space. See [8] or [9]. The Stratonovich form of (|47p 
is given by 

d s d t z st = d s d t w st - \d s z st dt, s,t^ 0. 
Fix x € M rf and consider for each t ^ the Stratonovich stochastic differential equation 

d s x st = Xi(x st )d s z l st + X (x st )ds, s ^ 0, 
with initial value Xot = This can be written in Ito form as 

(48) d s x st = Xi{x st )d s z l st + X (x st )ds, s ^ 0, 

where Xq = Xq + ^ X^i=i VXj.Xj. Consider also, for each t ^ 0, the stochastic differential 
equation 

d s U st = VXi(x st )U st d s zi t + VX (x st )U st ds, s > 0, 
with initial value Uot = I, and its Ito form 

(49) d s U st = VXi(x st )U st d s zi t + VX {x st )U st ds, s > 0. 

Proposition 5.1. There exist (two-parameter) semimartingales (z s t : s,t^ 0), (x s t : s,t ^ 
0) and {U s t ■ s,t ^ 0) snc/i £aat (z s t : s, t ^ 0) satisfies |^7[ ) and, /or aZZ t ^ 0, (x s t : s ^ 0) 
and (t/gt : s ^ 0) satisfy |^| ) and J^ffi ), mi/i £ae boundary conditions given above. Moreover, 
almost surely, U s t is invertible for all s,t ^ 0. 

Proof. We seek to apply Theorem 14.21 There are three minor obstacles: firstly to deal with 
the ds and dt differentials appearing in the equations, secondly, to show that the domain of 
the solutions is the whole of (M + ) 2 and, thirdly, to deal with the fact that the coefficients 
in (|49|) do not have the required boundedness of derivatives. 
Let us introduce a further equation 

d s d t z° st = 0, 

with boundary conditions z® = s and zot = t for all s, t ^ 0. We then replace dt and ds 
in (|4T|) and (08]), respectively, by dtZ% and d s z® t . When we obtain a solution, it will follow 
that Zg t = s + 1, so dtz® t = dt and d s z° st = ds, as required. 

In order to show that T> = (]R + ) 2 , it will suffice to show that the companion processes 
u s t and v st associated with the equations 

d s d t z° st = 0, d s d t z st = d s d t w st - \d s z st d t z Q st , 
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according to equations (jSJ) and (|1U|) . along with their inverses, remain bounded on compacts 
in s and t. We leave this to the reader. 

Finally, choose for each M £ N a smooth and compactly supported function tpM ° n 
R d <g) (R d )*, such that tp M (U) = U whenever \U\ < M. We can apply Theorem H21 to the 
system (|17|) . ([15]) . together with the modified equation 

Define 

£> M = : ^ M for all s' ^ s,t' ^ t}. 

By local uniqueness, we can define consistently U on P = Um'Dm by C/ S t = for 
€ £>a/. By some straightforward estimation using the one-parameter equations (|49l) . 
we obtain, for all T < oo and all p£ [1, oo), a constant C < oo such that 

sup E(\u st - u sltl \n {{sms , t , )m ) < c(| s - s ?/ 2 + \t - t'f 2 ). 

Then, by [9j Theorem 3.2.1], almost surely, U is bounded uniformly on T> n [0,T] 2 . Hence 
P = (]R + ) 2 , and we have obtained the desired semimartingale U . The invertibility of U can 
be proved by applying the same argument to the usual equation for the inverse. □ 

By the Stratonovich chain rule, 

d s d t x st = VXi(x st )d s z l st d t x s t + VX {x st )dsd t x st + Xi(x st )d s d t zl t . 

Now 

d s d t U s t = VXi{x st )d s z l st d t U s t + VX (x st )dsd t U st 

+ (^XiixstjdtXs^Ustdszlt + {V 2 X (x st )d t x st )U st ds + VXi(x st )U st d s d t zi t , 

so 

d t U st d s d t z l st = \d s d t Ustd s dtw l st = \VXi(x st )U st dsdt 

and 

d s (u^d t u st ) = u; t l {v^iixMzi^st + v 2 x (x st )dsd t x st + vx^d^} u st . 

Define also a two-parameter, Revalued, semimartingale (y s t ■ s,t ^ 0) by 

dtVst = U~ t l d t x sU y s0 = 0. 

Then 

d s d t y s t = U~ t l Xi{x st )d s d t z l st . 

Note that 

dty s td s dtz l st = d t y s td s dtw l st = \d s d t y st d s d t w % st = ^U^X^Xg^dsdt. 

So 

d s {d t y st <g> dtyst) = d s d t y st <8> d t y st + d t y st ® d s d t y st = U^X^Xst) <g> U'^X^Xs^dsdt. 
Note also that 

dsiU^Xiixst)) = U^[X u X 3 ]{x st )d s zi t + U-^Xo^ds. 

So 

d a (U* 1 X i (x llt ))d B d t 4 t = U* 1 [X i) X j ](x at )d s 4 t (d a d t wU - \d s z\ t dt) = 0. 
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Moreover 

dtiu^Xiix^dsdtzU = dtiu^Xiix^dsdtwit = o. 

Hence, we have 

d s dtyst = U'^X^Xs^dsdtzlt = U^ 1 X i (x st )(d s d t wl t - \d s z\ t dt). 

We compute 

dsiU^dtUstdtyst) 

= u; t l {^Xiix^dszU + v 2 x (x st )ds} d t x st ® d t x st + u^vx^xs^x^x^dsdt. 

Define 

R st = - U^XiiXrJdrZit, C st = U^XiiXr^^U^Xiix^dr. 

Jo Jo 
Our calculations show that the (J- S t ■ t ^ 0)-semimartingale (y s t : t ^ 0) has finite-variation 
part (y s t : t ^ 0) and quadratic variation given by 

dtyst = \Rstdt, dtyst ® 9 t y st = C st c2i. 

Moreover 

d t x st = U st d t y st + ld t U st d t y s t, 
so (x s t : t ^ 0) has finite-variation part (x s t : t ^ 0) and quadratic variation given by 

d t x st = \L st dt, d t x st <g> d t x st = T st dt, 

where 



r s 

Let = U xtR.lt + 



t f C/rtHV 2 ^^)^^ + V 2 X (x rt )5r}r rt 

JO 

+ U st [ U^-VXiix^Xiix^dr 
Jo 



and where T st = U st C st U* t . 

Note that both (T s t : t ^ 0) and (L s j : t ^ 0) are stationary processes and that, by 
standard one-parameter estimates, T s q and L s q have finite moments of all orders. By 
Ito's formula, for any C 2 function /, setting f st = f(x st ), the process (f s t : t ^ 0) is an 
(^" s t : i ^ 0)-semimartingale with finite-variation part (f s t : t ^ 0) and quadratic variation 
given by 

dtfst = \ (L st Vif{x st ) + rjjViVj/(zrt)) dt, a t / st a t / st = Vi/Cx^rjjVi/c^d*. 



In particular, if = / s j — / s q — /st) then (rn s t : t ^ 0) is a (true) martingale. Hence, for 
/,g £ C^(M d ), we obtain the integration- by-parts formula 

■ ' 1 
E[Vj/(x s0 )r l / Vj5(x s0 )] = 1 ™^ IE [{/(^at) - /(^so)}{5(a:st) - g(x s o)}] 

= -2lim^E[f(x s0 ){g(x st ) - g(x s0 )}] = -¥ l [f(x s0 ){Li V i g(x s0 )+F%V i V j g(x s0 )}}. 

An obvious limit argument allows us to deduce the following simple formula, corresponding 
to the case g(x) = xK For all / £ C 2 (M. d ) and for j = 1, . . . , d, we have 

E[V^(x s0 )lt ] = -E[/(x s0 )4 ]. 
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The general formula can then be recovered by replacing / by fVjg and summing over j. 

The basic observation underlying this formula is that the distributions of (zo,zt) and 
(z t ,zo) are identical, and hence that the same is true for (x S Q,x st ) and (x st ,x s0 ), when 
(x s t : s ^ 0) is obtained by solving a stochastic differential equation driven by (z st : s ^ 0), 
with initial condition independent of t. In fact a stronger notion of reversibility is true. 
The distributions of (z su :s)0,«e [0, i]) and (z Si t- u :s)0,ue [^M]) are identical, and 
hence the same is true for (x su : s ^ 0, u G [0,t]) and {x s ^- u '■ s ^ 0,u G [0, t]). This may 
be combined with the fact that the Stratonovich integral is invariant under time-reversal 
to see that 



E 



{f(x st ) - f(x s0 )} [ U-Jdux 
Jo 



-2E 



fix 



s0 



f 

Jo 



U RU d u x s 



From this identity, by a similar argument, we obtain the following alternative integration- 
by-parts formula. For all / £ C^(M d ), we have 

E[Vf(x s0 )U s0 C s0 } = -E[f(x s0 )R s0 ]. 

This formula is the variant discovered by Bismut, which is closely related to the Clark- 
Haussmann formula. 
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